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Existing limits on the non-radiative decay of one neutrino to another plus a massless particle (e.g., 
a singlet Majoron) are very weak. The best limits on the lifetime to mass ratio come from solar 
neutrino observations, and are r/m > 10 -4 s/eV for the relevant mass eigenstate(s). For lifetimes 
even several orders of magnitude longer, high-energy neutrinos from distant astrophysical sources 
would decay. This would strongly alter the flavor ratios from the 4> Ve : <(>„ : <j> Vr = 1:1:1 expected 
from oscillations alone, and should be readily visible in the near future in detectors such as IceCube. 


PACS numbers: 13.35.Hb, 14.60.Pq, 95.85.Ry 

Neutrinos from astrophysical sources are expected to 
arise dominantly from the decays of pions and their 
muon daughters, which results in initial flavor ratios 
(f) Ve : <j) v : (f > Vt of nearly 1:2:0. The fluxes of each mass 
eigenstate are given by <f> Vi =J2 a < /'t° urce |^ r ai| 2 , where U a i 
are elements of the neutrino mixing matrix. For three ac¬ 
tive neutrino species (as we assume throughout) there is 
now strong evidence to suggest that v tJ and v T are maxi¬ 
mally mixed and U e 3 ~ 0. The consequent Vu—v T symme¬ 
try means that in the mass eigenstate basis the neutrinos 
are produced in the ratios <fr Vl : (j) V2 : 0„ 3 = 1 : 1 : 1, in¬ 
dependent of the solar mixing angle. Oscillations do not 
change these proportions, but only the relative phases 
between mass eigenstates, which will be lost. An inco¬ 
herent mixture in the ratios 1 : 1 : 1 in the mass basis 
implies an equal mixture in any basis ( UIIA 4 = /), and 
in particular the flavor basis in which the neutrinos are 
detected (Ij. In this Letter we show that neutrino decay 
could alter the measured flavor ratios from the expected 
1 : 1 : 1 in a strong and distinctive fashion. 

We restrict our attention to two body decays 

Vi —> vj + X and v^ —> Fj + X, (1) 

where v, are neutrino mass eigenstates and X denotes 
a very light or massless particle, e.g. a Majoron. Vi¬ 
able Majoron models which feature large neutrino decay 
rates have been discussed in Ref. @]. We do not consider 
either radiative two-body decay modes (which are con¬ 
strained by photon appearance searches to have very long 
lifetimes j3}) or three-body decays of the form v —> vvv 
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(which are strongly constrained Q by bounds on anoma¬ 
lous Zvv couplings 0 )- In contrast, the limits on the 
decay modes considered here are very weak. Beacom 
and Bell have shown that the strongest reliable limit is 
r/m > 10~ 4 s/eV, set by the solar neutrino data |6]. 
This limit is based primarily on the non-distortion of the 
Super-Kamiokande spectrum Q, and takes into account 
the potentially competing distortions caused by oscilla¬ 
tions (see also Ref. |3]) as well as the appearance of active 
daughter neutrinos. It is very likely that the SN 1987A 
data place no limit at all on these neutrino decay modes, 
since decay of the lightest mass eigenstate is kinemati¬ 
cally forbidden, and even a reasonable v\ flux alone can 
account for the data 0, 9’j. 

The strongest lifetime limit is thus too weak to elimi¬ 
nate the possibility of astrophysical neutrino decay by a 
factor of about 10 7 x (L /100 Mpc) x (10 T eV/E) @. A 
few previous papers have considered aspects of the decay 
of high-energy astrophysical neutrinos. It has been noted 
that the disappearance of all states except v\ would pre¬ 
pare a beam that could in principle be used to measure 
elements of the neutrino mixing matrix, namely the ra¬ 
tios Uj, : : U 2 , 0. The possibility of measuring 

neutrino lifetimes over long baselines was mentioned in 
Ref. mt , and some predictions for decay in four-neutrino 
models were given in Ref. E2 We will show that the 
particular values and small uncertainties on the neutrino 
mixing parameters allow for the first time very distinctive 
signatures of the effects of neutrino decay on the detected 
flavor ratios. The expected increase in neutrino lifetime 
sensitivity (and corresponding anomalous neutrino cou¬ 
plings) by several orders of magnitude makes for a very 
interesting test of physics beyond the Standard Model; 
a discovery would mean physics much more exotic than 
neutrino mass and mixing alone. We will show that neu¬ 
trino decay cannot be mimicked by either different neu¬ 
trino flavor ratios at the source or other non-standard 
neutrino interactions. 




A characteristic feature of decay is its strong energy 
dependence: exp (— L/t\ & \,) = exp (—Lm/Er), where r 
is the rest-frame lifetime. However, we will assume that 
decays are always complete, i.e., that these exponential 
factors vanish. This is reasonable because there is a min¬ 
imum L/E value set by the shortest distances (typically 
hundreds of Mpc) and the maximum energies that will 
be visible in a given detector (the spectra considered 
are steeply falling). The assumption of complete decay 
means we do not have to consider the distance and inten¬ 
sity distributions of sources. We assume an isotropic dif¬ 
fuse flux of high-energy astrophysical neutrinos, and can 
thus neglect the angular deflection of d aug hter neutrinos 
from the trajectories of their parents ) 1 .‘ll |. It is uncer¬ 
tain if astrophysical sources produce the same numbers of 
neutrinos and antineutrinos. Though the detectors can¬ 
not distinguish neutrinos from antineutrinos, their cross 
sections are different, and this could cause confusion in 
the deduced flavor ratios. However, the antineutrino- 
neutrino cross section ratio is 0.7 at 10 TeV, and rapidly 
approaches unity at higher energies. 

Disappearance only. — We first assume that there 
are no detectable decay products, that is, the neutrinos 
simply disappear. This limit is interesting for decay to 
‘invisible’ daughters, such as a sterile neutrino, and also 
for decay to active daughters if the source spectrum falls 
sufficiently steeply with energy. In the latter case, the 
flux of daughters of degraded energy may make a negligi¬ 
ble contribution to the total flux at a given energy. Since 
coherence will be lost we have 

fo a {E) = J2K UICe ( E )\ U 0i\ 2 \ U ^\ 2 e~ L/Ti{E) (2) 

i(3 

4> S vp TCe (E)\U@i\ 2 \Uai\ 2 , (3) 

i(stable) ,(3 


where the (j> Va are the fluxes of v a , U a i are elements of the 
neutrino mixing matrix and are the neutrino lifetimes 
in the laboratory frame. Eq. m corresponds to the case 
where decay is complete by the time the neutrinos reach 
Earth, so only the stable states are included in the sum. 

The simplest case (and the most generic expectation) 
is a normal hierarchy in which both and vi decay, 
leaving only the lightest stable eigenstate v \. In this case 
the flavor ratio is U pl : t/ 2 -, : U%i H 3 ■ Thus if U e 3 = 0 
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FIG. 1: The effect of the presently unknown U e 3 on the 
<\>v e l<\> ratio. We have fixed 9q = 30° and 0 atm = 45°. Al¬ 
though varying these angles affects the flux ratios to a similar 
extent as U e 3 , they will be precisely measured in the near 
future. In all cases, the three scenarios are very distinct. 


V 3 is the lightest and hence stable state, and so 

foe : fo„ ■■ for = U* ■■ U% :U r 3 = 0 : 1 : 1. (5) 

If U e 3 = 0 and 9 atm — 45°, each mass eigenstate 
has equal v t j and v T components. Therefore, decay 
cannot break the equality between the </>„ and (f> VT 
fluxes and thus the 0 Ue : <j) v ratio contains all the use¬ 
ful information. The variation of the 0 Ve : fo ratio 
with non-zero U P ?. (up to the maximum allowed value, 
|C / e3 | 2 < 0.03 Q) is shown in Fig. 0 In the no-decay 
case, the variation from 1 : 1 : 1 is negligibly small. 
While the relative effect can be larger if neutrino de¬ 
cay occurs, the three cases shown are always quite dis¬ 
tinct. In addition, the ratio of the and v T components 
can also change, e.g., Eq. 0 could be as extreme as 
Uei '■ U?n '■ = 3.5 : 1 : 0.3. Hereafter, we set U e 3 = 0. 

Appearance of daughter neutrinos. — If neutrino 
masses are quasi-degenerate, the daughter neutrino car¬ 
ries nearly the full energy of the parent. An interesting 
and convenient feature of this case is that we can treat 
the effects of the daughters without making any assump¬ 
tions about the source spectra. Including daughters of 
full energy, we have 


foe '■ fo„ '■ for = c °s 2 ■ \ s in2 6*0 : \ sin 2 0 Q ~ 6 : 1 : 1 , 

(4) 

where 9@ is the solar neutrino mixing angle, which we 
have set to 30°. Note that this is an extreme deviation 
of the flavor ratio from that in the absence of decays. It 
is difficult to imagine other mechanisms that would lead 
to such a high ratio of u e to Here and throughout we 
concentrate on the flavor ratios, since the original source 
fluxes are unknown. In the case of an inverted hierarchy, 


fo{E) ^^ ource (£)|f/^| 2 |C/ QI | 2 

i/3 

+ Y. ( t )S r C ^E)\U 0J \ 2 \U ai \ 2 B^i ( 6 ) 

ij/3 

where B is a branching fraction and stable and unstable 
states are denoted henceforth by i and j respectively. 

If instead the neutrino mass spectrum is hierarchical, 
the daughter neutrinos will be degraded in energy with 
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respect to the parent, so that 


h *( E ) l ^| 2 |^| 2 (?) 

i(3 

/ oo 

dE' W E , E J2 4> s °r e {E')\U0j\ 2 \U m \ 2 B 3 ^, 
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where E is the daughter and E' is the parent energy. The 
normalized energy spectrum of the daughter is given by 


We'e 


1 dT{E',E) 
T {E') dE 


( 8 ) 


If the neutrinos are Majorana particles, daughters of 
both helicities will be detectable (as neutrinos or antineu¬ 
trinos), whereas if they are Dirac particles, daughters of 
one lrelicity will be sterile and hence undetectable. In 
the rest frame of the parent neutrino, the angular dis¬ 
tributions for decays which conserve and flip helicity are 
proportional to cos 2 (6* / 2) and sin 2 (#*/2) respectively, 
where 9* is the angle of the daughter neutrino with re¬ 
spect to the (lab frame) momentum of the parent. In the 
limit TOdaughter < ^parent, the corresponding energy dis¬ 
tributions in the lab frame are E/E' 2 and (E r — E)/E' 2 . 

In the case of Majorana neutrinos, we may drop the 
distinction between neutrino and antineutrino daughters 
and sum over helicities. Assuming the source spectrum 
to be a simple power law, E~ a , we find 


^ a (E) Ec; rce (^)i^i 2 i^r 

i(3 

+ ^EC irce (^i [/ «i 2 i^i 2 ^- (9) 

ij/3 

This is identical to the expression in Eq. © except for 
the overall factor of 1/a in front of the second term. For 
Dirac neutrinos we detect only the daughters that con¬ 
serve helicity, the effect of which is only to change the 
numerical coefficient of the second sum in Eq. © . Thus, 
although the flavor ratio will differ from the cases above, 
it is still independent of energy—i.e., decay does not in¬ 
troduce a spectral distortion of the power law. We stress 
that we have assumed a simple but reasonable power law 
spectrum E~ a ; a broken power law spectrum, e.g., would 
lead to a more complicated energy dependence. 

Uniqueness of decay signatures. — Depending on 
which of the mass eigenstates are unstable, the decay 
branching ratios, and the hierarchy of the neutrino mass 
eigenstates, quite different ratios result. For the normal 
hierarchy, some possibilities are shown in Table I. 

The most natural possibility with unstable neutrinos is 
that the heaviest two mass eigenstates both completely 
decay. The resulting flavor ratio is just that of the light¬ 
est mass eigenstate, independent of energy and whether 
daughters are detected or not. For normal and inverted 
hierarchies we obtained 6:1:1 and 0:1:1 respec¬ 
tively. Interestingly, both cases have extreme <p v< , : 


TABLE I: Flavor ratios for various decay scenarios. 


Unstable 

Daughters 

Branchings 

<^e ■ <t>v u ■ <t>VT 

V2, 

anything 

irrelevant 

6:1:1 

V3 

sterile 

irrelevant 

2:1:1 

V3 

full energy 
degraded (a = 2) 

= 1 

1.4 : 1 : 1 

1.6 : 1 : 1 

V3 

full energy 
degraded (a = 2) 

L>3_,1 = 1 

2.8 : 1 : 1 

2.4 : 1 : 1 

v 3 

anything 

B 3 ^i = 0.5 

B 3 ->2 = 0.5 

2:1:1 


ratios, which provides a very useful diagnostic. Assum¬ 
ing no new physics besides decay, a ratio greater than 1 
suggests the normal hierarchy, while a ratio smaller than 
1 suggests an inverted hierarchy. In the case that decays 
are not complete these trends still hold, even though the 
limits of Eqs. 0© would not be reached. The case of in¬ 
complete decay might be identified by measuring different 
flux ratios in different energy ranges. It is interesting to 
note that complete decay cannot reproduce 1:1:1. One 
of the mass eigenstates does have a flavor ratio similar 
to 1 : 1 : 1, but it is the heavier of the two solar states 
and cannot be the lightest, stable state. (A possible but 
unnatural exception occurs if only this state decays). 

An important issue is how unique decay signatures 
would be. Are there other scenarios (either non-standard 
astrophysics or neutrino properties) that would give sim¬ 
ilar ratios? There exist astrophysical neutrino produc¬ 
tion models with different initial flavor ratios, such as 
0:1:0 [HI, f° r which the detected flavor ratios (in the 
absence of decay) would be about 0.5 : 1 : 1. However, 
since the mixing angles 0 q and 0 atm are both large, and 
since the neutrinos are produced and detected in flavor 
states, no initial flavor ratio can result in a measured 
4 >v e '■ 4>v ratio anything like that of our two main cases, 
6:1:1 and 0:1:1. 

In terms of non-standard particle physics, decay is 
unique in the sense that it is “one-way”, unlike, say, oscil¬ 
lations or magnetic moment transitions. Since the initial 
flux ratio in the mass basis is 1 : 1 : 1, magnetic mo¬ 
ment transitions between (Majorana) mass eigenstates 
cannot alter this ratio, due to the symmetry between 
i —> j and j —► i transitions. On the other hand, if neu¬ 
trinos have Dirac masses, magnetic moment transitions 
(both diagonal and off-diagonal) turn active neutrinos 
into sterile states, so the same symmetry is not present. 
However, the process will not be complete in the same 
way as decay- it will average out at 1/2, so there is no 
way we could be left with a only single mass eigenstate. 

Experimental Detectability. — Deviations of the 
flavor ratios from 1:1:1 due to possible decays are 
so extreme that they should be readily identifiable flfij . 
Upcoming high energy neutrino experiments, such as Ice- 
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Cube a will not have perfect abilities to separately 
measure the neutrino flux in each flavor. However, the 
quantities we need are closely related to the observables, 
in particular in the limit of v )J —v T symmetry (0 atm = 45° 
and U e 3 = 0), in which all mass eigenstates contain equal 
fractions of and In that limit, the fluxes for and 
v T are always in the ratio 1 : 1 , with or without decay. 
This is useful since the v T flux is the hardest to measure. 

Detectors such as IceCube will be able to directly mea¬ 
sure the Vfj, flux by long-ranging muons which leave tracks 
through the detector. The charged-current interactions 
of u e produce electromagnetic showers. However, these 
may be hard to distinguish from hadronic showers caused 
by all flavors through their neutral-current interactions, 
or from the charged-current interactions of v T (an initial 
hadronic shower followed by either an electromagnetic or 
hadronic shower from the tau lepton decay) [13. We 
thus consider our only experimental information to be 
the number of muon tracks and the number of showers. 

The relative number of shower events to track events 
can be related to the most interesting quantity for testing 
decay scenarios, i.e., the v e to v fl ratio. The precision 
of the upcoming experiments should be good enough to 
test such extreme flavor ratios produced by decays. If 
electromagnetic and hadronic showers can be separated, 
then the precision will be even better. 

Comparing, for example, the standard flavor ratios of 
1 : 1 : 1 to the possible 6 : 1:1 generated by decay, 
the more numerous electron neutrino flux will result in a 
substantial increase in the number of showers compared 
to the number of muon events. The details of this ob¬ 
servation depends on the range of muons generated in 
or around the detector and the ratio of charged to neu¬ 
trino current cross sections. This measurement will be 
limited by the energy resolution of the detector and the 
ability to reduce the atmospheric neutrino background. 
The atmospheric background drops rapidly with energy 
and should be negligibly small above the PeV scale. 

Discussion and Conclusions. — We have presented 
our results above in terms of the ratios of fluxes in each 
neutrino flavor. These ratios are energy-independent be¬ 
cause we have assumed that the ratios at production 
are energy-independent, that all oscillations are aver¬ 
aged out, and that all possible decays are complete. The 
first two assumptions are rather generic, and the third 
is a reasonable simplifying assumption. In the standard 
scenario with only oscillations, the final flux ratios are 
4>v e '■ ■ <j>v T = 1:1:1. In the cases with decay, we 

have shown rather different possible flux ratios, for ex¬ 
ample 6 : 1 : 1 in the normal hierarchy and 0 : 1 : 1 in the 
inverted hierarchy. These deviations from 1:1:1 are so 
extreme that they should be readily measurable. 

These clear and striking predictions for the effects of 
neutrino decay on the measured flavor ratios depend 
strongly on recent progress in measuring neutrino mix¬ 
ing parameters. In particular, it is very significant that 


0q ~ 30° is well below the maximal 45°, for which 
Eq. Q would instead be a much less dramatic 2:1:1. 
In addition, 0q < 45° means that 5m\ 2 > 0 and hence 
that V 2 (with flavor ratios 0.7 : 1 : 1) can never be the 
lightest mass eigenstate. Maximal 0 atm and very small 
U e 3 also make the predictions clearer. The hierarchy of 
^3 relative to the two solar states is unknown, but in 
either case neutrino decay will be stringently tested by 
upcoming measurements of astrophysical neutrinos. 
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